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ABSTRACT 

A  "radiosity  correction"  algorithm  to  improve  the  quality  of  sidescan  sonar  mosaics  is 
reported.  The  algorithm  removes  consistent  range-dependent  variations  from  sidescan 
sonar  imagery  caused  by  failure  of  the  TVG  (time-varying  gain)  to  completely 
compensate  for  variations  due  to  transducer  beam  profiles  and  sediment  backscatter 
strength.  The  correction  is  appropriate  to  high-frequency  sidescan  sonars  with  constant 
gain  characteristics,  such  as  the  Klein  5000.  The  algorithm  appears  robust  in  situations 
where  the  bottom  is  relatively  flat  and  the  towfish  altitude  does  not  change  much. 
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Executive  Summary 

Sidescan  sonar  mosaics  are  a  geo-registered  (map-like)  form  of  seabed  imagery  with 
particular  utility  in  mine  hunting,  route  surveillance  and  hydrography.  They  can 
potentially  provide  high-quality  information  about  sediment  variations  and  object  and 
obstacle  locations.  This  report  describes  a  simple  "radiosity  correction"  algorithm  to 
improve  sidescan  sonar  mosaics  that  are  subject  to  consistent  range-dependent 
intensity  variations. 

The  algorithm  removes  consistent  range-dependent  variations  from  sidescan  sonar 
imagery  caused  by  failure  of  the  TVG  (time-varying  gain)  to  completely  compensate 
for  variations  due  to  transducer  beam  profiles  and  sediment  backscatter  strength.  The 
correction  is  appropriate  to  high-frequency  sidescan  sonars  with  constant  gain 
characteristics,  such  as  the  Klein  5000. 

The  algorithm  has  proven  effective  for  many  high-resolution  sidescan  sonar  collected 
with  a  Klein  5000  in  Sydney  Harbour  and  other  areas.  The  algorithm  appears  to  be 
very  robust  in  situations  where  the  bottom  is  relatively  flat  and  the  towfish  altitude 
does  not  change  much,  a  large  proportion  of  most  datasets.  The  algorithm  imposes 
minimal  processing  overheads  on  modem  personal  computers. 

The  algorithm  does  not  work  effectively  where  the  bottom  is  strongly  sloping.  It  can 
also  induce  shading  variations  that  are  not  directly  associated  with  sediment  type 
when  the  altitude  of  the  towfish  changes.  The  algorithm  is  not  recommended  for 
sediment  classification;  however,  sediment  boundaries  are  easier  to  see  in  mosaics 
prepared  using  the  algorithm. 
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1.  Introduction 

Geo-registered  or  "mosaiced"  sidescan  sonar  data  provides  a  map-like  view  of  the 
seabed  with  high  potential  value  for  applications  such  as  mine  countermeasures. 
However,  the  quality  of  a  mosaic  is  very  sensitive  to  defects  in  the  data  it  incorporates. 
One  of  the  defects  commonly  encountered  is  sometimes  referred  to  as  a  radiosity 
effect.  This  is  due  to  the  vertical  beam  pattern  of  the  transducers,  which  induce  a 
consistent  range-dependent  variation  in  the  signal  level. 

This  report  describes  a  simple  "radiosity  correction"  algorithm  that  has  proved 
effective  in  reducing  range-dependent  variations  in  imagery  from  the  Klein  bUUU 
sidescan  sonar.  The  Klein  5000  operates  at  455kHz.  It  is  a  short-range,  high-resoluhon 
device  with  a  relatively  high  signal  to  noise  ratio  across  most  of  the  swath,  a  factor  that 
is  important  to  the  success  of  the  algorithm.  It  also  has  fixed,  operator-independent 
gain  characteristics.  This  limits  the  number  of  different  correction  algorithms  that  must 

be  developed. 

The  radiosity  correction  algorithm  has  been  widely  applied  and  appears  to  be  quite 
robust.  It  is  not  greatly  affected  by  the  towing  altitude  of  the  towfish,  the  bottom  type 
or  the  water  depth.  However,  the  algorithm  has  not  been  tested  with  data  from  lower- 
frequency  sidescan  sonars  and  its  applicable  frequency  range  has  not  been  established. 


Figure  1  Sidescan  mosaic  before  (left)  and  after  (right)  correction  for  radiosity  variations 
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Figure  1  shows  an  example  of  the  effect  of  the  radiosity  correction.  The  bottom  has 
areas  of  silt,  shell-bed  and  rock.  Sediment  variations  are  suggested  by  differences  in 
shading  in  the  uncorrected  image,  but  the  cross-range  shading  variations  make  it 
difficult  to  identify  boundaries  between  them.  Sediment  boundaries  are  much  more 
readily  visible  in  the  corrected  image.  In  particular,  it  is  easy  to  see  the  exposed  shell- 
bed  that  marks  the  vessel  track  going  up  to  a  wharf  at  the  top-left  of  the  image. 

Mechanisms  responsible  for  radiosity  effects  are  discussed  in  Section  2  of  this  report, 
and  some  indication  is  given  as  to  where  the  technique  is  likely  to  succeed. 

The  implementation  and  processing  demands  of  the  algorithm  are  discussed  in  Section 
3,  along  with  improvements  that  could  be  made  in  future. 


2.  Radiosity  variations 


2.1  Intensity  variation  with  range 

The  total  electrical  response  of  the  transducers  can  be  described  in  terms  of  electrical 
power  as: 

PT  =t]Pa+Ne  (2.1) 

Here,  Pr  is  the  total  electrical  power  induced  in  the  system,  77  is  the  acoustic  efficiency 
of  the  transducers  and  amplifiers,  and  Ne  is  the  electrical  noise  power  added  by  the 
system.  After  amplification,  the  total  power  is 

Ps  =  g(t)PA+\0UFnog{t)NE  (2.2) 


Here,  g(t)  is  the  amplifier  gain  and  NF  is  the  noise  factor  of  the  amplifier  system.  Note 
that  g(t)  is  a  "time-varying  gain"  with  a  value  that  is  adjusted  in  an  attempt  to  maintain 
the  average  output  power  approximately  constant. 

The  acoustic  power  is  the  sum  of  "signal"  power  and  "noise"  power.  For  the  purposes 
of  bottom  characterisation,  the  signal  in  this  case  is  the  part  of  the  acoustic 
reverberation  from  the  seabed  that  follows  a  direct  path  to  and  from  the  transducers. 
All  other  acoustic  energy  reaching  the  transducers  is  "noise",  including  surface  and 
multipath  reverberation  and  environmental  noise.  Hence, 


P  -  ps  +  pN 

1  A  1  A  '  J  A 


(2.3) 


The  acoustic  signal  power  is  a  function  of  the  bottom  characteristics,  the  geometry  and 
the  beam  pattern  via 
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P  =  P 

rA  r0 


exp(-2ar)  2 


b2(o)sh(r) 


(2.3) 


Here,  the  time-independent  quantities  are  Po,  the  power  output  of  the  transducer  when 
it  emits  a  pulse;  a,  the  attenuation  rate  of  sound  in  water;  b,  the  transducer  beam 
pattern;  and  s>„  the  backscattering  strength  of  the  bottom.  The  others  are  time- 
dependent:  r=r(t),  the  one-way  range  (also  known  as  the  "slant  range");  O=0(t)  is  the 
start  angle  of  the  ray  that  will  reach  the  bottom  at  range  r,  relative  to  the  transducer 
surface  normal  -  the  "beam  angle";  and  y=y(t)  is  the  "grazing  angle"  that  the  ray  makes 
where  it  contacts  the  bottom. 

The  acoustic  noise  power  is  a  sum  over  contributions  from  the  surface,  ambient  sources 
and  multiple  bottom/ surface  reflected  energy  ("multipaths").  Each  of  these  terms  has  a 
form  with  similar  complexity  to  (2.3). 

Clearly,  the  power  recorded  by  the  system  is  a  complicated  function  of  the  signal  and 
multiple  noise  sources,  which  are  in  turn  complicated  functions  of  the  topography, 
environmental  characteristics  and  device  characteristics.  When  the  time-varying  gain  g 
of  the  system  is  imperfectly  matched  to  the  environment  and  the  system,  the  radiosity 
effects  seen  in  Figure  1  arise. 

A  correction  to  radiosity  variations  will  only  be  feasible  if  most  of  the  complicating 
effects  are  negligible.  The  regularity  of  the  effect  in  Figure  1  suggests  that  this  is  so. 

2.2  A  simplified  model 

If  we  assume  a  flat  sea  bottom  where  the  direct  reverberation  from  the  bottom  (the 
signal)  dominates  other  acoustic  inputs  to  the  transducers,  and  a  low-noise  sonar 
design,  then  the  total  electrical  power  recorded  by  the  system  will  be  approximated  by 

P  =vg«)P„  exP  (-2gr)^(eg(y)  (2.4) 

r 

If  we  further  assume  isovelocity  propagation  and  a  transducer  tilt  or  depression  angle 
of  8,  then  we  have 

0=Si"H(r)“5 

•  -1  f  Q  \ 

y  =  sin  — 


(2.5) 

(2.6) 
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Here,  the  altitude  of  the  transducers  is  a.  The  symbols  referred  to  in  equations  (2.3)  and 
following  are  depicted  in  Figure  2,  which  shows  a  transducer  at  altitude  a,  with  its 
surface  normal  depressed  by  an  angle  S.  A  point  on  the  bottom  at  slant  range  r  has  an 
angle  6  relative  to  the  transducer  face  and  sound  waves  strike  the  point  with  grazing 
angle  y. 


Figure  2:  Geometry  of  the  sidescan  swath,  looking  along-track 


2.2.1  Geometrical  dependencies 

The  rate  of  change  of  both  angles  is  strongly  range-dependent,  since 

dO  _  dy  _  - 1/  r  (2 

dr  dr  yjr2/a2  -1 

The  rate  of  change  is  infinite  at  the  nadir  point  where  a=r,  but  falls  rapidly  as  r  becomes 
much  larger  than  a.  Given  that  sidescan  sonars  are  deliberately  towed  close  to  the 
bottom,  this  is  the  case  over  much  of  the  swath. 

Figure  3  shows  the  variation  of  the  beam  angle  with  range  for  a  few  typical  towing 
heights.  Ranges  out  to  150m  are  displayed,  but  a  Klein  5000  might  typically  be  used 
with  a  75m  range  setting.  It  would  also  normally  be  deployed  less  than  20m  from  the 
bottom. 

Apart  from  the  rapid  change  of  beam  angle  near  nadir,  the  curves  in  Figure  3  show  that 
it  is  insensitive  to  range  once  the  range  has  exceeded  a  few  multiples  of  the  towfish 
altitude.  The  form  of  the  curves  also  becomes  similar,  and  they  approach  each  other  in 
value.  Finally,  if  the  towfish  is  close  to  the  bottom,  the  beam  angle  is  also  quite 
insensitive  to  the  altitude. 
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Figure  3  Variation  of  beam  angle  0  with  slant  range  r  from  the  sonar  for  some  typical  towfish 
altitudes.  The  transducer  depression  angle  is  10°  downward. 


Equations  (2.5)  and  (2.6)  imply  that  the  beam  angle  9  and  grazing  angle  y  have  the 
same  variation  with  range.  Grazing  angle  is  thus  insensitive  to  range  and  towfish 
altitude  when  the  towfish  is  run  close  to  the  bottom.  Hence,  at  long  ranges  from  the 
sonar,  a  single  radiosity  correction  may  be  valid  for  a  range  of  altitudes  if  the  beam 
pattern  and  backscatter  strength  vary  only  slowly  with  angle.  This  is  unlikely  to  be  the 
case  closer  to  the  sonar. 

Figure  3  also  implies  that  variations  due  to  altitude  will  result  in  a  constant  shading 
offset  at  longer  ranges,  since  the  gradients  of  the  curves  are  quite  similar. 

2.2.2  Backscatter  variation 

Backscatter  models  valid  at  455kHz  have  not  appeared  in  the  literature  to  date.  An 
empirical  model  published  by  the  University  of  Texas  [McKinney  and  Anderson  1964] 
is  claimed  to  be  valid  for  the  frequencies  from  20- 290kHz.  A  semi-empirical  model 
published  by  the  Applied  Physics  Laboratory  at  the  University  of  Washington  [APL 
1984, 1989, 1994]  is  claimed  to  be  valid  in  the  more  restricted  range  10-100kHz. 
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Figure  4  APL  model  predictions  of  bottom  backscatter  strength  for  generic  sediment  types  at 
4551<Hz. 

In  order  to  proceed,  we  must  choose  one  of  the  models  and  use  it  well  outside  of  its 
domain  of  validity.  Given  that  the  APL  model  is  physically  based,  we  assume  that  it 
will  predict  relative  behaviour  reasonably,  even  if  the  absolute  values  for  the 
backscatter  strength  are  incorrect. 

Figure  4  shows  bottom  backscatter  strength  as  a  function  of  grazing  angle  for  a  range 
of  "generic"  sediment  types.  It  is  notable  that  at  this  high  frequency,  the  curves  do  not 
cross,  excepting  at  the  highest  grazing  angles.  All  of  the  curves  excepting  sandy  gravel 
have  similar  forms  for  grazing  angles  less  than  about  70°,  implying  that  the 
backscattered  intensity  at  a  given  angle  should  be  a  reasonably  good  indicator  of 
sediment  grain  size. 
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Figure  5  Bottom  backscatter  strength  for  455kHz  sound  on  a  generic  medium  sand  bottom  as  a 
function  of  range  along  the  bottom  and  towfish  altitude. 


Figure  5  shows  bottom  backscatter  strength  as  a  function  of  horizontal  range  along  the 
bottom  for  a  generic  "medium  sand"  bottom,  at  455kHz.  The  predicted  backscatter 
strength  varies  strongly  as  a  function  of  both  horizontal  range  and  towfish  altitude. 
The  appearance  of  Figure  4  suggests  that  we  can  remove  the  variation  with  range 
almost  independently  of  the  sediment  type.  However,  the  variation  with  altitude 
would  be  more  difficult  to  remove  and  would  tend  to  make  shallow  areas  appear 
acoustically  less  reflective.  This  means  that  a  mosaic  based  on  a  range-dependent 
radiosity  correction  would  not  be  useful  for  sediment  classification  unless  the  towfish 
were  flown  at  a  constant  altitude  from  the  bottom.  A  corrected  mosaic  may  however 
still  be  useful  for  the  location  of  sediment  boundaries. 
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3.  Radiosity  correction  algorithm 

3.1  Processing  order 

The  first  decision  that  must  be  made  when  applying  a  radiosity  correction  is  when  to 
apply  it,  in  particular,  should  the  algorithm  be  applied  before  or  after  the  sidescan  data 
is  projected  from  its  original  dependence  on  slant  range  (equivalently,  time)  to 
horizontal  range  across  the  bottom,  so-called  "slant-range  corrected"  data. 

Figure  4  shows  that  the  intensity  variation  in  a  sidescan  sonar  image  is  likely  to  be 
more  rapid  in  the  region  of  the  nadir  (high  grazing  angles)  than  elsewhere.  The 
proportion  of  the  image  affected  by  this  rapid  variation  is  small  (typically  about  5%) 
when  the  data  is  projected  as  a  function  of  the  slant  range  and  somewhat  larger  when  it 
is  projected  as  a  function  of  horizontal  range.  For  this  reason,  it  was  decided  to  apply 
the  radiosity  correction  to  the  data  before  slant-range  correction. 

In  the  radiosity  correction  algorithm,  all  calculations  start  at  the  sample  corresponding 
to  the  first  bottom  return  -  they  are  "bottom-tracking"  and  rely  on  an  accurate 
determination  of  the  location  of  the  first  bottom  return.  Samples  corresponding  to 
returns  from  the  water  column  are  ignored. 

3.2  The  filtration  process 

No  single-stage  filtration  process  investigated  by  the  author  appeared  to  be  sufficient 
to  remove  range  variations.  Invariably,  some  variation  remained  in  the  near-nadir 
region.  Hence,  a  two-stage  process  was  adopted.  First,  variations  away  from  the  nadir 
region  are  removed  with  a  fixed  correction.  Then,  an  adaptive  correction  is  applied  to 
remove  the  remaining  variation  near  nadir. 

3.2.1  First  stage 

Considering  only  one  side  of  the  sidescan  data,  let  the  time  series  of  Ns  samples  for 
ping  n  be 

s(n,i),  i  =  0,l...Ns-\  (3-1) 

Let  the  first  bottom  return  for  ping  n  be  sample  i  =  b(n).  We  assume  that  a  process 
exists  for  estimating  b(n)  accurately. 

Let  A/min  be  the  smallest  number  of  'on-land'  samples  per  ping  encountered  during  the 
averaging  process,  which  is  given  by 

Nmm  =  mm(Ns  - b(n)),  n  =  \..NP  (3.2) 
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Clearly,  if  a  range-dependent  variation  is  consistent  from  ping  to  ping,  then  it  should 
become  evident  after  averaging  sufficiently  many  pings  along  track.  However,  the 
variation  will  follow  the  position  of  the  first  bottom  return.  Thus,  we  take  a  "bottom¬ 
tracking"  average  of  the  signal  along  the  track,  via 

a(i)  =  ^-^s(n,i-b(n)),i  =  0.  .Nmm  (3.3) 

np 


Nr  is  the  total  number  of  pings  averaged  along  track.  A  value  of  Nr  =  200  has  been 
found  to  be  adequate  in  this  case. 


Note  that  the  average  specified  by  Nr  must  be  long  enough  to  remove  local  variations, 
but  not  long  enough  to  go  from  one  major  area  of  seabed  sediment  type  to  another. 


After  the  first  average,  a(i)  remains  a  relatively  unsmooth  function.  When  applied 
directly,  it  causes  streaks  to  appear  along-track  in  the  corrected  image,  corresponding 
to  spikes  in  a(i).  To  remove  the  spikes,  a  rolling  average  is  applied  across-track  to 
smooth  the  function  further,  via 


i  '/(o 

a{i)  =  — - - 7T  7  Z,aU) 

qi) )  -  pi1) + 1 W 


Here, 

p(i)  -  max(0,  i  -  /) 
q(i)  =  min(yVmin  -1,/  +  /) 
l  =  Umj50 


(3.4) 


(3.5) 


Finally,  the  average  signal  level  for  the  entire  area  is  evaluated.  Let  the  across-track 
average  of  the  along-track  averages  be  A,  given  by 


A=- 


1 


N, 


min 


(3.6) 


min 


The  first-stage  correction  is  given  by 

C(i)  =  A /  max(l,  a  (/))  (3.7) 

The  correction  is  applied  multiplicatively  via 

sc  ( n ,  i )  -C(i-  b(n))  s(n,  i)  (3.8) 


9 


DSTO-TN-0354 


Note  again  that  the  correction  is  tied  to  the  first  bottom  return.  Additive  corrections  do 
not  work,  since  the  beam  pattern  variation  itself  has  a  multiplicative  effect  on  the  signa 
level.  Note  also  that  the  correction  has  no  effect  on  the  signal  in  the  water  column. 

Figures  6  and  7  show  the  effects  of  the  first  stage  correction  on  adjacent  sections  of  data 
taken  from  one  of  the  runs  shown  in  Figure  1.  Figure  6  shows  along-ttack  aver ages 
of  the  port  and  starboard  signal  levels  resulting  from  equation  (3  3).  Note  that  0  slant 
range  fs  the  slant  range  at  the  first  bottom  return.  There  is  clearly  a  large,  systema  , 
range-dependent  effect  on  both  channels.  The  across-track  smoothing  effect  of  equation 
(3.4)  on  the  port  and  starboard  averages  is  shown  by  the  cyan  and  green  lines 

respectively. 

Figure  7  shows  the  effect  of  applying  the  smoothed  correction  shown  in  Figure  6  to  an 
adjacent  section  of  data,  via  equations  (3.7)  and  (3.8).  Most  of  the  systematic  variation 
with  range  has  been  removed.  Only  a  rapid  variation  in  the  first  3  metres  remains. 


range,  m 


Figure  6  Sonar  envelope  voltages,  averaged  along-track,  before  correction 
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range,  m 


Figure  7  Sonar  envelope  voltages ,  averaged  along-track,  after  first-stage  correction 

3.2.2  Second  stage  -  near-nadir  correction 

After  the  initial  correction  has  been  applied,  the  image  becomes  more  uniform,  except 
in  the  near-nadir  region  where  very  rapid  variations  remain.  A  more  local,  adaptive 
correction  is  appropriate  to  this  section  of  the  image. 

The  near-nadir  correction  is  based  on  the  requirements  that,  at  the  point  where  the  port 
and  starboard  images  meet,  they  should  have  the  same  value  and  derivative.  In  this 
case,  both  'value'  and  'derivative'  refer  to  time  averages  over  scales  appropriate  to 
tonal  variations.  It  is  therefore  assumed  that,  close  to  the  nadir,  the  average  variation  of 
the  signal  level  away  from  a  straight  line  is  a  measure  of  the  gain  mismatch  at  that 

point. 

In  this  case,  the  'straight  line'  is  taken  to  be  the  line  between  two  'anchor'  points 
arbitrarily  located  1/10*  of  the  sonar  range  setting  distant  from  the  first  bottom  return. 
The  value  at  each  anchor  point  is  estimated  using  a  smoothing  average  along-track,  via 

v*  =  —((£- 1)  +  5*  (n,  b(n)  +  A))  (3-9) 

L 
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Here,  L  is  the  length  of  the  filter  and  A  is  the  number  of  samples  corresponding  to 
l/10th  of  the  range  setting.  The  '+'  and  symbols  refer  to  port  and  starboard 
respectively,  numbering  samples  outward  from  the  towfish.  An  appropriate  value  of  L 
has  been  experimentally  determined  to  be  30,  corresponding  to  movement  of  order  6m 
along- track. 

In  order  to  retain  the  detail  in  the  near-nadir  region,  the  same  kind  of  along-track, 
smoothing  average  is  applied  to  the  data  between  the  two  anchor  points.  Let  the  value 
of  each  along-track  average  be  vn(i),  given  by 

v*(i)  =  y((^-l)vw_|(i)  +  s±(«,W  +  0)  '  =  0..A-1  (3.10) 

Lj 


The  near-nadir  correction  is  given  by 


Ak,+(v,7-/0- 


/v,7(0  ' =  0..A 


AT,  =  «  +v;)/2 


This  is  applied  multiplicatively,  via 
Sc  («,  0  =  c;7  (i  -  b(n ))  sc  (w,  i) 


(3.11) 


(3.12) 


The  final,  corrected  signal  has  been  derived  from  the  initially  corrected  signal  by  a 
multiplicative  correction  that  does  not  act  on  the  water  column.  Note  that  the  value  of 
the  smoothing  length  L  must  be  adjusted  so  that  it  is  equivalent  to  several  multiples  of 
the  largest  expected  details  in  the  near-nadir  region,  but  should  otherwise  be  kept  as 
small  as  possible. 

Figure  8  shows  along-track  signal  level  averages  for  approximately  the  same  section  of 
data  as  Figure  7,  after  application  of  the  second-stage  correction  given  by  equations 
(3.11)  and  (3.12).  Longer  ranges  are  unaffected  by  the  correction,  but  the  rapid  change 
in  the  level  in  the  first  few  metres  of  range  has  been  removed.  The  port  and  starboard 
signal  levels  now  approximately  agree  in  both  value  and  derivative. 

The  left  and  right  halves  of  Figure  1  show  a  mosaic  of  the  same  data,  before  and  after 
both  stages  of  the  correction  algorithm  have  been  applied.  Clearly,  the  algorithm  has 
been  successful  in  this  case  in  reducing  the  tonal  variation  at  the  centre  of  each  line, 
and  across  the  mosaic.  The  four  lines  in  the  mosaic  may  still  be  distinguished,  but  the 
underlying  sediment  variations  (which  have  been  confirmed  by  Roxann  echo-sounder 
analysis  and  grab-sampling  [Hamilton,  2000])  are  clearly  visible. 
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r^jrs^wr^«i=^jrs:: 

clipped  region  due  to  the  emission  of  each  ping,  and  me  clear  a 
are  the  returns  from  the  water  column. 

isstsslsfli 

highest  at  intermediate  ranges.  The  variation 

ss =swsir- “ 

horizontal  plane  representing  true  range  . 

bottom  returns,  but  does  not  have  the  tonal  fall-off  at  long  ranges,  or  the  peak  at  mid 
ranges. 


13 


DSTO-TN-0354 


After  stage  2  of  the  correction,  the  final  image  appears  like  the  right  third  of  Figure  9, 
which  has  even  tonal  variation  at  all  ranges. 


Figure  9  Sidescan  sonar  raster  data,  uncorrected  (left),  after  stage  1  (centre)  and  after  stage  2 
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4.  Limitations 

The  radiosity  correction  in  this  report  has  been  tested  with  sets  of  data  from  various 
versions  of  the  Klein  5500  sidescan  sonar.  It  has  been  successful  where  data  was 
collected  over  a  relatively  flat  seabed. 

The  primary  limitation  of  the  correction  is  that  it  cannot  normalise  tonal  variations  over 
sloping  ground,  particularly  when  runs  are  not  parallel  to  each  other.  The  backscatter 
coefficient  in  equation  (2.4)  is  dependent  on  the  grazing  angle  that  the  sediment  patch 
makes  relative  to  the  sonar,  and  this  will  generally  be  different  when  the  sonar  views  a 
sloping  area  from  two  different  aspects. 

A  second  limitation  is  the  depedence  of  the  shading  on  towfish  altitude.  The  radiosity 
correction  does  not  attempt  to  correct  for  the  reduction  in  signal  level  when  the  towfish 
is  flown  closer  to  the  bottom.  Hence,  mosaics  made  with  this  technique  cannot  be  used 
for  sediment  classification  unless  the  towfish  altitude  has  been  held  constant. 

A  third  limitation  is  due  to  the  pointwise  nature  of  the  algorithm,  which  was  adopted 
to  minimise  processing  time.  Were  a  faster  processor  available,  it  would  be  better  to 
express  the  correction  as  a  polynomial  function  of  either  range  or  grazing  angle.  This 
would  avoid  problems  in  very  shallow  water,  where  the  number  Nm in  of  samples 
defined  in  equation  (3.3)  can  be  somewhat  larger  than  it  is  when  the  towfish  is  flown 
higher.  In  effect,  a  shallow-water  correction  can  be  used  for  deeper  water,  but  not  vice- 
versa. 

Figure  10  shows  a  mosaic  incorporating  the  flat  area  shown  in  Figure  1  at  the  left  side 
of  the  image,  but  also  incorporating  stongly  sloping  regions.  The  mosaic  is  aligned 
with  north  up  in  the  usual  way.  The  areas  to  the  east,  west  and  southeast  of  the  central 
peninsula  are  quite  flat,  but  a  deep  hole  is  located  in  the  south  and  southwest.  In  this 
sloping  region,  clear  tonal  variations  can  be  seen  between  runs,  sufficient  to  obscure 
features  of  interest.  Elsewhere,  tonal  variations  are  due  to  sediment  changes  from  shell- 
beds  (dark)  to  mud  (light). 

Note  that  some  of  the  tonal  variation  in  the  shallow  water  to  the  northeast  may  be  due 
to  a  reduction  in  towfish  altitude,  but  the  sediment  is  silty  and  hence  has  a  low 
backscatter  strength. 

The  technique  has  not  been  trialled  with  100kHz  data,  which  may  exhibit  a  less 
uniform  variation  of  backscatter  coefficient  with  sediment  type  than  is  observed  at 
455kHz.  If  backscatter  strength  curves  for  different  sediments  cross  at  low  or  medium 
grazing  angles,  then  the  correction  will  fail  on  flat  ground. 
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in  ^P^re^i0n 


5.  Conclusion 

This  report  has  demonstrated  a  simple  technique  for  removal  of  : 

artifacts  from  mosaics  prepared  using  Klein! SOW '  The  correction 

so  that  they  are  more  usehr.  for  sediment 

classification  purposes. 

before  slant-range  correction. 

with  a  constant  towfish  altitude. 
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